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© Deformation measuring method and device using cross-corre.atlon function between speck.e 
patterns. 



© In a deformation measuring method and device 
in which an object is irradiated with a laser beam 
before and after deformation of the object to obtain 
speckle patterns, the speckle patterns thus obtained 
are photoelectrically converted into electrical signals, 
and the cross-correlation function between the 
speckle patterns is .calculated using the electrical 
signals to obtain displacement of the speckle pattern 
on the basis of the shift of position of the extreme 
value of the mutual-correlation function and to deter- 
mine the amount of deformation of the object from 
the displacement of the speckle pattern, (1) the 
reference speckle pattern data for calculation of the 
cross-correlation function is renewed when the ex- 
treme value of the cross-correlation function is lower 
in level than a predetermined value or when the 
position of the extreme value is out of a predeter- 
mined range, or (2) the reference speckle pattern 
data is renewed when the extreme value of the 
cross-correlation function is lower in level than a 
predetermined value and a calculation range for cal- 
Oculating the cross-correlation function is shifted when 
Q/he position of the extreme value is out of a pre- 
UJ determined range. 
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formation measurement method according to the 
invention will be described with reference to FIG. 2. 

As shown in FIG. 2, a laser beam 13 outputted 
by a laser source 12 is applied to a measuring 
region 0 of an object 10. through a magnifying lens 
14 if necessary, and the resultant speckle pattern is 
observed on an observation plane 30. In this case, 
it is assumed that the coordinate axes on the 
object's surface are represented by x. y and z; the 
distance OS between the divergence point of the 
laser beam 1 3 and the measuring region O, Ls (OS 
= Ls)- components representing the direction, from 
the divergence point to the center of the region 0. 
Isx, isy and isz: the distance between the ob- 
ject's surface and the observation plane 30, Lo; 
Components representing the direction from the 
center of the region O to an observation point P, 
tx ly and Iz; and the components of translation, 
rotation and distortion of the object 10 at the region 
irradiated by the laser beam 13. (ax. ay, az). (Ox. 
lly, flz) and («xx, «yx, «yy). respectively. 

Speckle patterns are formed at the observation 
point P before and after deformation oMhe object. 
The cross-correlation function C( x . y ) between 
the intensity distributions li(x.y) and l 2 (x,y) of the 
speckle patterns thus formed is as follows: 
C( 5^ y- } = <|, ( X ,y) x l 2 (x + x Y + y )> (1) 
where <> means a set average. 

Calculating equation (1 ), it can be understood 
that C( x~, y~) has a maximum value with x = Ax 
and y~= Ay, where Ax and Ay are represented by 
the following equations (2) and (3) and physically 
correspond to the amount of displacement of the 
speckle pattern caused by the deformation of the 
object: 

Ax =-ax[(Lo/Ls)(tsx 2 - 1) + lx 2 -1] 

-ay[(Lo/Ls) tsx Isy + tx ly] 

-azf(lo/Ls) tsx isz + tx iz] - Lo[-Qz(?.sy + ly) 

-ay(lsz '+ iz) + £ xx(lsx + tx) + «xy(lsy + ty)l 

(2) 

Ay =-ax[(Lo/Ls)(lsy Isx + ly lx] 
-ay[(Lo/Ls)(isy 2 - D + ly 2 -1] 
-azKLo/Ls) Isy isz + lx Iz] - Lo[-Qz(tsx + lx) 
-fix(isz + tz) + e yy(lsy + iy) + «xy(lsx + lx)] 
(3) 

Accordingly, when the amount of displacement 
of the speckle pattern (Ax and Ay) are observed 
with one-dimensional image sensor disposed at the 
observation plane 30. the. output waveform of the 
image sensor changes as shown in FIG. 3A after 
and before deformation of the object, the auto- 
correlation waveform is as shown in FIG 3(B), and 
the cross-correlation waveform is as shown in FIG. 

3(C) 

In measuring the deformation of an object with 
the above-described device, uniike the convention- 
al method in which the reference data (speckle 
pattern) is fixed or renewed every time to obtain 
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the cross-correlation function, the reference data 
(speckle pattern) is renewed in this invention only 
when the extreme value of the cross-correlation 
function is lower than a predetermined value or the 
position of the extreme value comes out of a pre- 
determined range. Thus, when the speckle pattern 
is displaced greatly with the deformation of the 
object, the reference pattern is renewed and there- 
fore the range of measurement is not limited. In 
contrast, in the case where the deformation of the 
object is small and therefore the position of the 
extreme value is within the predetermined range, 
according to the method of the invention the refer- 
ence speckle pattern is fixed, so that the displace- 
ment at a low speed can be detected- unlike the 
conventional method in which the reference data is 
renewed every time, accordingly, the accuracy of 
measurement can be increased while the range of 
displacement measurement is increased with no 
errors. 

FIG. 4 shows the arrangement of one example 
of a deformation measuring device to which the 
deformation measuring method according to the 
invention is applied. The device, as shown in FIG. 
4, comprises: a laser source 12 for applying a laser 
beam 13 to the surface of an object 10 to be 
measured to form a speckle pattern. In order to 
subject the speckle pattern to photo-electric con- 
version, two Photosensitive element arrays 20 and 
24 are arranged with an angle therebetween. Each 
of the photosensitive element arrays has a strip or 
comb structure in which each of the photosensitive 
elements (-, 22 n .„ 22 ni 22„*,. -) constituting a 
photosensitive element array has a rectangular 
form of a large ratio of a long shide to a short side 
(for example, a short side of 1 3 urn and a long side 
of 2.5 mm) are arranged at intervals of 25 urn in a 
strip form as shown in FIG. 5. 

The speckle patterns, after subjected to pho- 
toelectric conversion by the photosensitive element 
arrays 20 and 24, are applied to correlators 40 and 
42 according to the invention, in each of which the 
change in position of the extreme value of the 
cross-correlation function before and after the 
speckle pattern displacement is.detected. 

The correlators 40 and 42 are identical in con- 
struction, and therefore only the correlator 40 will 
be described with reference to FIG. 6. 

The correlator 40 as shown in FIG. 6. com- 
prises: a memory 40A for temporarily storing the 
preceding speckle pattern (one-frame) data of the 
binary-coded (A/D-converted) electrical signal A 
which is obtained by the preceding scanning op- 
eration through the photosensitive element array 20 
(in the case of the correlator 42, through the pho- 
tosensitive element array 24): correlation IC 40B 
(for instance TDC1023 manufactured by TRW 
Company) for receiving a reference speckle pattern 
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data B from the memory 40A and a comparison 
speckle pattern data C from the photosensitive 
element array 20 or 24 to output the correlation 
value D between the reference and comparison 
data 8 and C; a gate 40C for renewing the refer- 
ence speckle pattern data B; a maximum value 
detecting circuit 40D for outputting the extreme 
value E of the correlation values D provided every 
shift clock pulse CL, and its timing output F; a 
counter 40F for receiving the timing output F and 
the shift clock pulse CL to output a signal G 
representing the position of the extreme value (in 
this case, maximum value): and a comparison cir- 
cuit 40F for comparing the extreme value E output- 
ted from the maximum value detecting circuit 40D 
with a predetermined threshold value, to output a 
reference pattern exchange signal H when the ex- 
treme value E is less than the threshold value. The 
correlator 40 further comprises: a comparison cir- 
cuit 40G for comparing the shift of the position G of 
the extreme value outputted from the counter 40E 
with a predetermined maximum threshold value for 
the maximum range, and outputting a reference 
pattern exchange signal I when the position of the 
extreme value exceeds the maximum threshold val- 
ue- a comparison circuit 40H for comparing the 
shift of the position G of the extreme value with a 
predetermined minimum threshold value (for the 
minimum range), and outputting a reference pattern 
exchange signal J when the position G of the 
extreme value is less than the minimum threshold 
value; an Or gate 40I for performing add operation 
between the reference pattern exchange signals H, 
I and J and outputting the resultant to the gate 40C 
and an FIFO (first-in first-out) memory 40J 
(described later): and the FIFO memory 40J for 
temporarily storing the the position G of the ex- 
treme value outputted from the counter 40E and 
the reference pattern exchange signal H + l + J and 
applying them to a CPU bus 40K. 

According to a flow chart shown in FIG. 7, 
when the extreme value is lower in level than the 
threshold value therefor or the position G of the 
extreme value is out of the predetermined range, 
the renewing operation of the reference pattern is 
automatically carried out. 

The shift of the position of the extreme values 
detected by the correlators 40 and 42 are applied 
to and processed by the computer 44, in which the 
amount of deformation of the object is calculated 
from the amount of displacement of the speckle 

pattern. ■ . „• 

The computer 44 applies an instruction to the 

stepping motor controller 46 to move along the x- 

axis the linear stage 48 on which the object 10 is 

placed so that the speckle pattern is displaced 

along the x-axis. and at the same time applies 

necessary timing signals to the correlators 40 and 



42. 

FIG. 8 is a graph showing the relation between 
the measured and actual data with the threshold 
value of the extreme value E set to 100%. 90%. 
s 80%. 70% and 60% in the case where the linear 
stage 48 is moved 6 mm along the x-axis. In the 
measurement using the threshold value of 100%. 
the reference speckle pattern is renewed every 
time This is identical to the conventional measure- 
io ment in which a speckle pattern of the just preced- 
ing frame is employed as the reference pattern. In 
the measurement using the threshold value of 
60% the reference pattern is not renewed because 
the threshold value is low. This is substantially 
,s identical to the conventional measurement in which 
the reference pattern is fixed. 

On the other hand, it is apparent from FIG. 8. 
in the case where the reference speckle pattern is 
renewed with the threshold values ranged from 
20 70% to 90%, the measured displacement of the 
linear stage which is obtained from the displace- 
ment of the speckle pattern corresponds substan- 
tially to the actual movement of the linear stage 48. 
In the first embodiment of the method for ob- 
25 taining the cross-correlation function as described 
above, when the deformation of the object is made 
in the direction intersecting to the observing- plane 
of the object, that is. the same position of the 
object is irradiated by the laser beam before and 
so after deformation, the speckle pattern is moved 
over the image sensor with small displacement. In 
this case, with the extreme value of the cross- 
correlation function less decreased, the position of 
the extreme value is moved. Accordingly, the num- 
35 ber of times at which the position of the extreme 
value comes out of the predetermined range is 
increased and therefore the number of exchanging 
the reference pattern of the cross-correlation pat- 
tern is greatly increased. An exchange of the refer- 
40 ence pattern causes an error of the half of the pitch 
interval of the photosensitive element in the image 
sensor at maximum, so that the total error cor- 
responding to the product of the pitch interval and 
the number of exchanging the reference pattern 

45 may be caused. 

In order to overcome this difficulty, the follow- 
ing second embodiment is employed: When, in 
obtaining the cross-correlation function, the ex- 
treme value of the cross-correlation function is 

so smaller than a predetermined value, the reference 
pattern of the mutual-correlation function is re- 
newed: and when the position of the extreme value 
comes out of the predetermined range, the refer- 
ence pattern of the cross-correlation function is not 

ss renewed, and instead thereof the range of calcula- 
tion for cross-correlation function is changed (or 
shifted). Hence, when the same part of the object 
is irradiated by the laser beam before and after 
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deformation of the object as in the case where the 
object is displaced out of the observing plane of 
the object, the frequency of exchanging the refer- 
ence pattern is decreased and therefore the mea- 
surement error can be decreased. Accordingly, in 
measurement of the displacement of an object 
when the speckle pattern is less changed, the error 
is minimized and the accuracy of measurement is 
increased. 

The deformation measuring device to which the 
second embodiment of this invention is applied has 
the same fundamental construction as that in the 
first embodiment shown in FIG. 4; however, it 
should be noted that the correlators of the first and 
second embodiments are greatly different from 
each other. 

The correlator 40 (or 42) of the second em- 
bodiment as shown in FIG. 9. comprises three shift 
registers (140A, 140B and 140C) for holding refer- 
ence pattern data B for the binary-coded electrical 
signal A of the speckle pattern provided by the 
photosensitive element array 20 or 24), for holding 
comparison data C for calculation of the correlation 
function with the reference pattern data B, and for 
holding the binary-coded electrical signal A pro- 
vided by the photosensitive element array 20 (or 
24) during calculation of the correlation function. 

The roles of three shift registers are dynam- 
ically swapped for one another by means of a 
sequencer 140D; that is, the shift registers are 
used alternately as a reference pattern. data shift 
register, a comparison data shift register, and a 

buffer shift register. 

The shift register selected for holding the refer- 
, ence pattern data is connected in ring mode by the 
sequencer 140D, so that the reference pattern data 
is held for the next correlation calculation. 

The sequencer 1400 selects the input of a 
multiplexer (MPX) 140E to apply one of the outputs 
of the three shift registers, as the reference pattern 
data B, to a correlation calculation circuit 140G; 
and further it selects the input of the other mul- 
tiplexer (MPX) 140F to apply another of the outputs 
of the three shift registers, as the comparison data 

C, to the circuit 140G. 

The correlation calculation circuit 140G is 
shown in detail in FIG. 10. The circuit 140G com- 
prises: a shift register 140H; an array 1401 of. for 
example, 16 EXCLUSIVE OR circuits; and an array 
140J of 16 binary counters, so that sixteen points 
of the correlation function can be simultaneously 
calculated. 

A maximum value detecting circuit 140K re- 
ceives the values of sixteen points of the correla- 
tion function which has been calculated by the 
correlation calculation circuit 140G. to detect the 
extreme value of the correlation function and its 
position E. and outputs them to comparison circuits 



140L, 140M and an FIFO (first-in first-out) memory 
140N. 

In the comparison circuit 140L the extreme 
value D provided by the maximum value detecting 
s circuit 140K is compared with a predetermined 
threshold value Dth. When the extreme value D is 
lower then the threshold value Dth, the comparison 
circuit 140L applies a reference pattern exchange 
signal F to the sequencer 140D and the FIFO 

w memory 140N. 

On the other hand, in the comparison circuit 
140M, the position E of the extreme value provided 
by the maximum value detecting circuit 140K is 
compared with a predetermined correlation calcula- 
75 tion range threshold value H. According to the 
result of the comparison, the comparison circuit 
140M applies a correlation calculation range shift 
signal I to the sequencer 140D and the FIFO mem- 
ory 140N to shift the correlation calculation range 
20 to the positive side or to the negative side. 

In response to the reference pattern exchange 
signal F. the sequencer 140D causes the roles of 
the shift registers 140A, 140B and 140C to be 
exchanged for one another, and changes the se- 
as lected input terminals of the multiplexers 140E and 
140F. In addition, in response to the correlation 
calculation range shift signal I, the sequencer 140D 
applies a shift clock pulse to the reference pattern 
data shift register to shift the correlation calculation 

30 range. 

As shown in a flow chart of FIG. 11, when the 
extreme value D becomes lower than the predeter- 
mined threshold value Dth, the reference pattern is 
automatically renewed and when the position E of 

35 the extreme value comes out of the predetermined 
range H, the correlation calculation range is auto- 
matically shifted. 

As shown in FIG. 4, the shifts in the position of 
the extreme value detected by the correlators 40 

40 and 42 are applied to and processed by the com- 
puter 44, and the amount of deformation of the 
object to be detected is calculated from the amount 
of displacement of the speckle pattern. 

The computer 44 applies an instruction to the 

45 stepping motor controller 46 to move the linear 
stage 48 along the x-axis on which the object 10 is 
placed, thereby to displace the speckle pattern 
along the x-axis, and at the same time applies 
timing signals to the correlators 40 and 42. 

so In the above-described first and second em- 

bodiments, the photosensitive element arrays 20 
and 24 comprising photosensitive elements are in 
the strip (comb) form and each photosensitive ele- 
ment has a rectangular form of a large ratio of a 

55 long side to a short side, that is, is in a strip of 
fancy paper, and further at least two of the pho- 
tosensitive element arrays are provided. As a re- 
sult, it is unnecessary to coincide the direction of 
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displacement of the speckle pattern with the direc- 
tion of the photosensitive element array, and the 
components in each of the directions can be mea- 

The configuration of the photosensitive ele- 
ments and the number of photosensitive element 
arrays are not limited only to those in the above- 
described embodiments. For instance, in the case 
where it is required to extract only the components 
of a speckle pattern along the x-axis, the number of 
photosensitive element arrays to be provided may 
be only one. In addition, the provision of more than 
two photosensitive element arrays will allow the 
measurement with high accuracy. 



Claims 

1. A deformation measuring method for detect- 
ing the amount of deformation of an object from 
displacement of a speckle pattern of the object, 
comprising the steps of: 

irradiating the object with a laser beam before and 
after deformation of the object to obtain speckle 

patterns; . . 

converting the speckle patterns into electrical sig- 
nals; 

calculating cross-correlation function between the 
speckle patterns from the electrical signals and 
obtaining the extreme value of the cross-correlation 
function; 

obtaining the amount of position shift, of the ex- 
treme value of the cross-correlation function to 
determine the amount of the deformation of the 
object, wherein said calculating step comprises the 
steps of: comparing a comparison speckle pattern 
data with a reference speckle pattern data for the 
cross-correlation function and renewing the refer- 
ence speckle pattern data when the extreme value 
of the cross-correlation function is lower in level 
than a predetermined value. 

2. A deformation measuring method as claimed 
in claim 1, wherein said renewed reference speckle 
pattern data comprises a speckle pattern data ob- 
tained by the just preceding scanning operation to 
the present scanning operation providing the com- 
parison speckle pattern data. 

3. A deformation measuring method as claimed 
in claim 1, wherein said calculating step further 
comprises the step of renewing the reference 
speckle pattern data when the position of the ex- 
treme value of the cross-correlation function comes 
out of a predetermined range. 

4. A deformation measuring method as claimed 
in claim 3, wherein said predetermined range has 
the minimum value corresponding to a half pitch 
interval of a photosensitive element array provided 
for performing said converting step. 



5. A deformation measuring method as claimed 
in claim 1, wherein said calculating step further 
comprises the step of shifting the calculation range 
for the cross-correlation function when the position 

5 of the extreme value of the cross-correlation func- 
tion is out of a predetermined range. 

6. A deformation measuring method as claimed 
in claim 5, wherein said shifting step comprises the 
step of shifting the calculation range at a positive 

70 side or a negative side thereof. 

7. A deformation measuring device for detect- 
ing the amount of deformation of an object from 
displacement of a speckle pattern of the object 
comprises: 

75 a light source for irradiating the object with a laser 
beam to obtain speckle patterns before and after 
deformation of the object; 

at least one photosensitive element array compris- 
ing plural photosensitive elements for converting 
20 the speckle patterns of the object into the cor- 
responding electrical signals; 

at least one correlator for calculating the cross- 
correlation function between the speckle patterns 
from the electrical signals and obtaining the posi- 
25 tion of the extreme value of the cross-correlation 
function, said correlator comprising memory means 
for temporarily storing the preceding speckle pat- 
tern data, correlation calculating means for receiv- 
ing a reference speckle pattern data from said 
30 memory and a comparison speckle pattern data 
from said photosensitive element array and output- 
ting the correlation values between the reference 
and 'comparison speckle pattern data, a maximum 
value detecting circuit for outputting the extreme 
35 value of the correlation values, comparison means 
for comparing the level of the extreme value and 
the position thereof with a predetermined value and 
a predetermined range respectively and outputting 
an instruction signal for renewing the reference 
40. speckle pattern data when the extreme value is 
less in level than the predetermined value or when 
the position of the extreme value is out of the 
predetermined range, and a gate for renewing the 
reference speckle pattern data to the preceding 
45 speckle pattern data in response to the instruction 
signal; 

a microcomputer for obtaining the amount of shift 
in position of the extreme value of the cross-cor- 
relation function and determining the amount of 
so deformation of the object from the amount of the 

shift . , . . 

8. A deformation measuring device as claimed 

in claim 7, wherein said photosensitive element 
array has a comb structure. 
55 9. A deformation measuring device for detect- 

ing the amount of deformation of an object from 
displacement of a speckle pattern of the object, 
comprises: 
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a light source for irradiating the object with a laser 
beam to obtain speckle patterns before and after 
deformation of the object: 

at least one photosensitive element array compris- 
ing plural photosensitive elements for converting 5 
the speckle patterns of the object into the cor- 
respondent electrical signals; 
at least one correlator for calculating the cross- 
correlation function between the speckle patterns 
from the electrical signals and obtaining the posi- w 
tion of the extreme value of the cross-correlation 
function, said correlator comprising three shift reg- 
isters for serving as a first shift register for holding 
the electrical signals as reference speckle pattern 
data, as a second shift register for the electrical 
signals as comparison speckle pattern data and a 
buffer shift register for holding the electrical signals 
during calculation of the cross-correlation function, 
said three shift registers being dynamically 
swapped in function for one another, a first mul- so 
tiplexer for receiving output signals from said shift 
registers and outputting one of the output signals 
as the reference speckle pattern data and a second 
multiplexer for receiving output signals from said 
shift registers and outputting one of the output as 
signals as the comparison speckle pattern data, a 
correlation calculation circuit for receiving the refer- 
ence and comparison speckle pattern data from 
said first and second multiplexers and calculating 
the cross-correlation function therebetween, a ao 
maximum value detecting circuit for obtaining the 
extreme value of the cross-correlation function and 
the position thereof, a first comparison circuit for 
comparing the extreme value with a predetermined 
value and outputting a first instruction signal for 
renewing the reference speckle pattern data when 
the extreme value is lower in level than the pre- 
determined value, a second comparison circuit for 
comparing the position of the extreme value with a 
predetermined range and outputting a second in- 
struction signal for shifting the cross-correlation 
range for calculation of the cross-correlation func- 
tion when the position of the extreme value is out 
of the predetermined range, and a sequencer for 
exchanging the functions of said three shift regis- 
ters for one another and changing the selected 
input terminals of said first and second multiplexers 
in response to said first and second instruction 
signals to thereby renew the reference speckle 
pattern data and shift the cross-correlation calcula- so 
tion range. 

10. A deformation measuring device as 
claimed in claim 9, wherein said photosensitive 
element array has a comb structure. 
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© Deformation measuring method and device using cross-correlation function between speckle 
patterns. 



© in a deformation measuring method and device 
in which an object is irradiated with a laser beam 
before and after deformation of the object to obtain 
speckle patterns, the speckle patterns thus obtained 
are photoelectrical^ converted into electrical signals, 
and the cross-correlation function between the 
speckle patterns is calculated using the electrical 
signals to obtain displacement of the speckle pattern 
on the basis of the shift of position of the extreme 
value of the mutual-correlation function and to deter- 
mine the amount of deformation of the object from 
the displacement of the speckle pattern, (1) the 
reference speckle pattern data for calculation of the 
P) cross-correlation function is renewed when the ex- 
treme value of the cross-correlation function is lower 
fgin level than a predetermined value or when the 
f- position of the extreme value is out of a predeter- 
mined range, or (2) the reference speckle pattern 
f^data is renewed when the extreme value of the 
^cross-correlation function is lower in level than a 
w predetermined value and a calculation range for cal- 
©culating the cross-correlation function is shifted when 
the position of the extreme value is out of a pre- 
qj determined range. 
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